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Abstract 

The condensation transition, leading to complete mutual synchronization in 
large populations of globally coupled chaotic Rossler oscillators, is investi- 
gated. Statistical properties of this transition and the cluster structure of 

partially condensed states are analyzed. 
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Appearance of condensates is an important property of quantum and classical systems. 
When a population of interacting Bose particles (such as, e.g., superfluid liquid He 4 ) is 
considered, a condensate is formed by a finite fraction of particles occupying the same 
quantum state Condensation of electron- hole pairs in metals is responsible for the 
phenomenon of superconductivity ||. Condensation phenomena in classical systems are 
closely related to the effect of mutual synchronization of oscillations. In large populations of 
globally coupled limit-cycle oscillators, one or several groups of oscillators following exactly 
the same phase trajectory can spontaneously appear fH-0. A similar behavior in populations 
of globally coupled logistic or circle maps has been found by K. Kaneko f|||. Mutual 
synchronization of several coupled chaotic oscillators has been studied (see, e.g., |IJ). To 
analyze properties of the condensation transition, which involves progressive formation of 
coherent oscillator groups, larger systems must be however considered. In this Letter we 
investigate condensation phenomena in a population of a thousand globally coupled chaotic 
Rossler oscillators. 

Not every form of global coupling between elements can lead to mutual synchronization. 
Synchronous chaotic oscillations are possible if 'vector' coupling (with identical coupling co- 
efficients for all oscillator variables) is employed [|TTJ. Asymptotic synchronization conditions 
in the limit of high coupling intensities have also been constructed for the case when cou- 



pling involves only one of the oscillator variables |T(J • However, only the linear stability of 
synchronous oscillations is proven in both of these cases and therefore even for a strong cou- 
pling the presence of other dynamical attractors, corresponding to persistent asynchronous 
oscillations, cannot be excluded ||12|| . 

In this Letter a different form of global coupling that yields more robust synchroniza- 
tion is introduced. We take an arbitrary dynamical system with m variables w(t) = 
{wi(t), w m (t)} whose evolution is determined by equations 

w = f(w), (1) 

where f(w) are some known functions. A globally coupled population consisting of iV such 



dynamical systems (i = 1, 2, N) is constructed by adding the coupling terms, so that the 
resulting evolution equations are 

wt = f(w») + £A(w - w<) + e'[f(w) - f (wO], (2) 

where w is the global average, w(i) = A^ -1 ^i=i w i(^); the positive parameters e and e' 
specify the intensity of global coupling, and A is a constant matrix. Note that the previously 
used forms of coupling correspond to a choice e' = 0; the 'vector' coupling is obtained for 
A = I where I is the identity matrix. 
When e' — 1 equations (@) reduce to 

Wj = eA(w — Wj) + f (w). (3) 

They describe linear motion of elements forming the population under action of a force which 
is collectively produced by all of them. The evolution of the deviations Swi(t) = Wj(t) — w(t) 
is then determined by an exact linear set of equations 5wj = — eA^Wj. Therefore, these 
deviations exponentially decrease with time, and global stability of synchronous oscillations 
is thus guaranteed, if all real eigenvalues of the matrix A and the real parts of all its complex 
eigenvalues are positive. 

Since globally stable synchronous oscilations are present in the system @ at e' = 1, 
they may also be expected in a neighborhood of this point, i.e. for an interval of coupling 
intensities Eq < e' < 1. As shown below, this interval can be so wide that even weak coupling 
leads to the synchronization. 

We consider a population of chaotic Rossler oscillators. Each elementary oscillator is 
described now by three variables w(t) = {x(t),y(t), z(t)}, and we have f x = —y — z, f y = 
x + ay, and f z — b — cz + xz, where a, b and c are fixed parameters. The collective dynamics 
of a globally coupled population of identical Rossler oscillators (i = 1, 2, N) is governed 
by equations 
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iji = Xi + ayi + K(y -yi), ( 4 ) 
ii = b - czi + XiZi + K(xz - XiZi), 

where x(t), y(t) and z(t) are global averages, such as x(t) = iV" 1 ^i(^), and K is the 
intensity of global coupling. This form of coupling is derived from (|2|) by taking e = e' = K 
and choosing a matrix A whose elements are zero except for A xy = A xz = —A yx = 1, 
A zz = c, and A yy = 1 — a. The system (|J) reduces to a set of independent identical Rossler 
oscillators ||13| for K = 0. 

The eigenvalues of the matrix A are Ai = c and A2,3 = \ (1 — a) ± — a ) 2 — 4. 

All three eigenvalues (or their real parts) are positive if a < 1. Therefore, for K = 1 all 
deviations from the global averages exponentially decrease with time and the states of all 
elements in the population asymptotically converge. We see that at K — 1 the oscillator 
population (|j) possesses a global attractor which corresponds to coherent motion of all 
elements of the system. This attractor coincides with the attractor of a single Rossler 
oscillator and is chaotic for the parameters a = 0.15, b = 0.4 and c = 8.5 used in our 
simulations. Thus, the population (||) of globally coupled chaotic Rossler oscillators should 
undergo condensation as the coupling intensity K is gradually increased towards K — 1. 
To investigate the condensation phenomena, we have performed numerical simulations of a 
population of iV = 1000 such oscillators under varying intensity K of global coupling. Below 
we first analyze statistical properties of the condensation transition and later consider how 
condensation is influenced by heterogeneities. 

To characterize condensation, distances = [(xj — Xj) 2 + (yi — yj) 2 + (zi — Zj) 2 } 1 / 2 between 
the states of all iV(iV — l)/2 possible pairs of elements have been computed. When 
condensation occurs, some of the elements have exactly the same coordinates and distances 
between the states of these elements are zero. The ratio of the number of pairs with zero 
distances to the total number of pairs can thus be chosen as the order parameter r of the 
condensation transition. In absence of a condensate, r = 0. On the other hand, r = 1 when 
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complete mutual synchronization of the whole population takes place. 

To prepare the initial condition, the system was first allowed to evolve up to t = 100 
without global coupling (K = 0) so that its elements get uniformly distributed over the 
Rossler attractor. Global coupling was then introduced, time was reset to zero, and the 
evolution of the coupled population was started. During a transient T, identical element 
pairs were appearing. The transient was over when the number of these pairs did not further 
increase with time. In our simulations, the transients were always shorter than T = 1000. 

Figure 1 shows the dependence of the order parameter r on the global coupling intensity 
K obtained by averaging over 20 realizations at t = 2000. Condensation begins at K c ps 0.017 
when a nonvanishing fraction of identical pairs first appears. In the interval from K c to 
K ps 0.06, nonmonotonous dependence of the order parameter on the coupling intensity is 
observed. At higher coupling intensities K, the fraction of identical pairs steadily grows until 
full condensation (r = 1) is established for K > K , with K ps 0.093. Under conditions of 
partial condensation, the number of identical pairs varies from one realization to another. 
Bars in Fig. 1 indicate statistical dispersion of the simulation data at a few selected values 
of the global coupling intensity. Note that, once a certain number of identical pairs has been 
formed at the end of the transient period, it remains later constant in any given simulation. 
Therefore, the dispersion reveals a degeneracy in the asymptotic properties of this system 
in the partially condensed state. 

The order parameter r characterizes the overall degree of condensation, but is not sen- 
sitive to the condensate structure. To determine the detailed structure of the population, 
distributions over pair distances between the states of all elements at a fixed time moment 
(t = 2000) have been analyzed. Histograms of such distributions at four different coupling 
intensities K are shown in Fig. 2. They are constructed by counting pairs with distances 
lying within subsequent intervals of width Ad = 0.25. The number of pairs inside each next 
interval, divided by the total number of pairs N(N — l)/2, yields the height of the respec- 
tive bar. In absence of global coupling, elements are relatively uniformly distributed over 
the single-oscillator attractor and a smooth distribution over distances d is thus observed 
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(Fig. 2a). When global coupling is introduced but remains below the critical point K c of 
the condensation transition, the distribution is modified (Fig. 2b). Now, pairs with small 
distances have already appeared, though identical pairs are still absent (cf. Fig. 1). Slightly 
above the critical point, a strongly nonuniform distribution with many narrow peaks is seen 
(Fig. 2c). The distributions found at larger intensities of global coupling (e.g., Fig. 2d) are 
formed by several distinct lines. When complete condensation has taken place (r = 1), the 
histogram has only one line located at d — (not shown). Distributions formed by several 
lines are characteristic for situations when the whole population breaks down into a number 
of coherent clusters. All elements in a cluster follow the same dynamical trajectory and, 
hence, distances between pairs of elements belonging to the same cluster are zero. When 
elements from two different clusters are chosen, their distances are all equal. If M such 
clusters are present, the distribution over pair distances would consist of M[M — l)/2 + 1 
individual lines. Hence, four lines seen in Fig. 2d correspond to a condensate with three 
clusters. In our simulations, pure three-cluster regimes were typically found starting from 
K 0.045. They were replaced by pure two-cluster regimes at K ~ 0.058. 

Hence, the condensate appearing above the transition point includes only a small fraction 
of elements of the population, so that the order parameter is small in this region. The 
condensed elements are distributed over a large number of small clusters. The rest of the 
population has asynchronous dynamics, though it is already essentially influenced by global 
coupling between the elements. When the coupling intensity is increased, the number of 
clusters decreases. However, the typical sizes of clusters become at the same time larger 
and therefore the condensate contains a greater fraction of the population. Starting from 
a certain intensity of global coupling, all elements belong to one of a few big clusters. As 
the coupling intensity is further increased, relative sizes of clusters vary and some of them 
disappear. This process eventually leads to complete condensation which can be viewed as 
a single-cluster regime. 

Fig. 3 shows the typical dependence of the order parameter r of the condensation tran- 
sition on the total population size N, where each point is obtained by averaging over 20 
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independent realizations and intensity of global coupling is kept constant (K = 0.08). We 
see that a strong size dependence is characteristic for relatively small populations (N < 500). 
Starting from a population size of about N = 1000, the dependence displays saturation, with 
the remaining small variations lying within the statistical dispersion range. Hence, the re- 
sults of our study of a population with 1000 oscillators might already be representative for 
the asymptotic behavior in the infinite-size limit iV — > oo. 

The complete condensation (r = 1) has been found within a wide interval of the global 
coupling intensity. Starting in this region with random initial conditions for all elements, 
we see that the elements soon form a compact cloud that rapidly shrinks with time. The 
characteristic radius R(t) of this cloud can be defined by 

N 

where 5xi, 5yi and 8zi are the deviations from the respective global averages. The radius 
R decreases with time, as the trajectories of all elements asymptotically converge to the 
same orbit (in our simulations the convergence was followed until the trajectories of all 
elements became identical up to the computer precision of about 10~ 6 ). An important 
quantitative property of this regime is the characteristic time which the system needs to 
reach the condensed state; the inverse of this time represents the condensation exponent 



1 N 

R 2 (t) = *f E N (*) + SvKt) + szm] , (5) 



(closely related to the transverse Lyapunov exponents discussed in Ref. flO|l ): 



7 = - lim - In — H 6 

The simulations show that 7 starts to grow from zero at K = Kq, has a maximum at 
K w 0.65 and then slowly decreases approaching the value 7 ~ 0.38 at K = 1. 

The robustness of the complete condensation suggests that it might persist in some form 
even in heterogeneous populations. To test this, we have carried out simulations where 
the elements forming the population were not identical. Heterogeneity was introduced by 
replacing the constant parameter c in Eqs. (|J) by random numbers q that were uniformly 
distributed inside the segment [c — a c ,c+ a c }. Alternatively, the same procedure was applied 



to the parameter a. Starting from random initial conditions, a compact cloud of elements 
was again formed in the heterogeneous case. However, the cloud did not shrink until it 
transformed into a single point. Instead, after an initial decrese, its radius R(t) fluctuated 
around a certain mean value. Fig. 4 shows the mean radius (R) as function of dispersions a c 
or a a for K = 0.2. We see that in a wide interval of heterogeneities the mean radius of the 
cloud is approximately linear proportional to the dispersion. Moreover, even for relatively 
large dispersions this radius is significantly smaller that the characteristic attractor diameter 
D, representing the maximal possible pair distance at K — (i.e. D 25, cf. Fig. 2a). 
Therefore, even though the states of all elements are no longer identical (i.e. the complete 
condensation is destroyed), their variations remain small. The elements form a little cloud 
whose center moves along a definite trajectory. 

This behavior can be understood by analyzing the effect of weak heterogeneities at the 
maximal strength of global coupling K — 1. If the parameter c is replaced by q = c + Aq 
with random variations Aq, the population dynamics is described by 

%i Hi ^"ii 

Vi = Xi + (a- l)yi + y, (7) 
— b z% I x z . 

Hence, motion of each element is governed by linear equations of motion in the presence of 
global driving forces y and xz , which are the same for all elements. Therefore, deviations 
Szi from the global average obey the equation 

5zi = —c5zi — AciZ — AciSzi + AciSzi. (8) 

If heterogeneities are weak, i.e. Aq = 8Q where is a small parameter, the last two terms in 
this equation have order 2 and can be neglected. Integrating the resulting linear equation, 
we obtain 

SzJt) = 0d / exp(-cr)^(t - r) dr. (9) 
J o 

Considering the behavior of other deviations 5xi and 5yi, one can similarly show that they 
are also proportional to Aq = 0Q. This means that the radius (R) of the population cloud 
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depends as (R) ~ 9 on the heterogeneity strength. Note also that in this case all deviations 
are scaling-identical, i.e., for example, Szi(t) = 9QQ(t) where the function Q(t) is the same 
for any element i. The same results are obtained when weak random variations of the 
parameter a are considered. 

Our simulations have revealed that the linear dependence of the mean radius (R) on the 
dispersions a c and a a , that has been analytically shown above only at K — 1, persists in 
a wide range of coupling intensities. Moreover, this dependence is found even for relatively 
large dispersions, where the parameter 6 would not be small. This behavior is apparently a 
consequence of a strong compression provided in the entire interval of complete condensation 
by global coupling. 

The condensation transition has been studied above for a particular kind of chaotic 
oscillators and for a certain form of global coupling. We have however checked that similar 
results are also obtained when a different ('vector') form of global coupling is used. This 
suggests that observed statistical properties may be typical for various large populations of 
globally coupled chaotic oscillators undergoing a condensation transition. 
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FIGURES 

FIG. 1. Order parameter r of the condensation transition as function of the coupling constant 
K. Averaging over 20 independent realizations is performed; bars indicate statistical dispersion of 
data at several selected points. 

FIG. 2. Normalized histograms of the distribution over pair distances d between elements in a 
population of N = 1000 globally coupled Rossler oscillators at different stages of the condensation 
transition: (a) K = 0; (b) K = 0.01; (c) K = 0.02; (d) K = 0.05. 

FIG. 3. Dependence of the order parameter r on the population size N at a fixed intensity of 
global coupling K = 0.08. Averaging over 20 independent realizations; bars indicate statistical 
dispersion of data. 

FIG. 4. Radius (R) of the cloud formed by a heterogeneous population as function of dispersions 
a a (empty circles) and a c (filled circles) at K = 0.2. 



11 



0.3 



0.2 



0.1 



a 



0.0 ""' —m***mmm 

0.3 

0.2 - 



0.1 



0.0 



Li. 



o 



10 



20 



d 




qp 
jo 
N 



CD 

Rp 
2 



o 

/ — - 

3 

f2 







10 



20 



d 



0.8 




10' 



10 







<R> 



10 



-1 







10 



-2 



10 



-3 




